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FLUTTER ANALYSIS OF TWO PARALLEL ELASTICALLY 
COUPLED FLAT PLATES 


Charles P. Shore 
Langley Research Center 

SUMMARY 

Flutter of two parallel elastically coupled flat plates is investigated analytically. 

A closed-form solution including both aerodynamic and structural damping is presented 
for flutter of flat orthotropic plates coupled by an elastic medium. Both plates are sim- 
ply supported along the side edges but are supported by deflectional, rotational, and tor- 
sional springs of arbitrary stiffness at the leading and trailing edges. Two-dimensional 
quasi-steady aerodynamics is utilized in the solution. Since the large number of vari- 
ables present in the problem precludes extensive parametric studies, results are pre- 
sented to indicate the basic flutter characteristics of coupled two-plate systems and to 
assess the validity of previously published modal solutions for similar problems. 

INTRODUCTION 

Flutter of two simply supported parallel isotropic plates connected by an elastic 
medium has been investigated previously in connection with micrometeroid bumpers 
which are attached to a primary structural skin by a light soft filler material. Although 
intended primarily as a protective device in space, such configurations must withstand 
severe aerodynamic loads during launch. In reference 1 a two-mode Galerkin solution 
to the coupled fourth-order differential equations describing motion of the two plates was 
obtained. The simple Ackeret expression for the lateral loading due to air pressure was 
used in the analysis. In reference 2 a two-mode flutter solution which used piston theory 
aerodynamics to incorporate aerodynamic damping effects was obtained for the coupled- 
plate problem. In addition a semiexact procedure and a finite-difference formulation 
were outlined; however, the only results presented were from the two-mode solution. 
Results presented in references 1 and 2 indicate that the elastic coupling medium can 
have significant effects on the flutter response of the two-plate system. In fact, depending 
on the magnitude of inplane loads acting on the plates and the spring stiffness of the elastic 
coupling medium, the two-plate system can become unstable at much lower dynamic pres- 
sures than either of the two plates alone. 

More recently several thermal protection systems which employ exposed surfaces 
coupled to primary load carrying structures by flexible insulation and/or strain-isolator 


systems have been proposed for space shuttle and hypersonic cruise vehicles. The flut- 
ter characteristics of reusable surface insulation (RSI) panels coupled to a primary struc- 
ture by a flexible strain isolator have been studied in reference 3. The analysis utilizes 
plate modes for the simply supported primary structure and free-free beam modes for 
the RSI panels. Dynamic response of the coupled system was obtained by application of 
the Rayleigh-Ritz technique. Results presented in reference 3 indicate that the presence 
of the RSI panels has a highly stabilizii^ effect on the flutter response of the primary 
structure. The main cause of the increased flutter stability of the primary structure 
was attributed to the segmented mode shapes of the RSI panels which govern the interac- 
tion of the aerodynamic forces and the deflections of the primary structure. Finally, the 
typical severely destabilizing influence of midplane compressive loads on single isotropic 
panels was not always observed for the coupled system, and the buckling value of the mid- 
plane loading for the primary structure was found to be increased significantly by the 
RSI/strain-isolator system. 

In many instances the exposed surfaces of typical thermal protection systems con- 
sist of highly orthotropic panels with arbitrary flexible supports at the boundaries. (See 
ref. 4.) Thus, the present report extends the previous solutions to a more general con- 
figuration consisting of two parallel orthotropic plates connected by an elastic medium 
and supported at the leading and trailing e<^es on springs of arbitrary stiffness. Piston 
theory aerodynamics is employed to obtain the lateral loading due to air pressure over 
the exposed plate, and taking the side edges of the plates to be simply supported allows a 
closed- form solution (i.e., nonmodal) to the governing partial differential equations. 

Results from the analysis are compared with results from references 1 to 3, and a gen- 
eral flutter boundary for two identical simply supported coupled plates is presented for a 
single value of the elastic spring stiffness of the connecting medium. Specific examples 
are also presented to indicate the applicability of the solution and to indicate some effects 
of structural damping and inplane loading on the flutter of elastically coupled parallel 
plates. Additionally, the influence on flutter of single plates with unequal support springs 
at the leading and trailing edges is briefly discussed. 

SYMBOLS 

Dimensional quantities are presented both in the International System of Units (SI) 
and in the U.S. Customary Units. Measurements were made in U.S. Customary Units. 

Aj deflection shape coefficients 

A parameter defined by equation (16) 

a plate length 
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D. 


^xy 

Di,D2 
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^x>^xy>^y 




Sk 


^x’^xy ’^y 
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i 

K 
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k^,k|-,kr 


parameter defined by equation (17) 
plate width 

free-stream speed of sound 
plate bending stiffness in x-direction 
plate bending stiffness in y-direction 
plate twisting stiffness 

plate stiffness coefficients defined by equation (7) 

Young ^s modulus 

complex deflectional, torsional, and rotational support spring coeffi- 
cients defined by equation (11) 

complex coefficients defined by equation (8) 

aerodynamic damping coefficient, ^ 

deflectional, torsional, and rotational support spring structural damping 
coefficients 

elastic coupling- medium structural damping coefficient 
plate structural damping coefficients 
moment of inertia 
imaginary unit, 

complex spring stiffness of elastic coupling medium 
spring stiffness of elastic coupling medium in compression 
deflectional, torsional, and rotational spring stiffnesses 
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Z(x,y,t) 


aerodynamic loading term 


M 

m 

Nx 

Ny 

n 

Pi 

q 

Rx 

RxjCr 

Ry 

s 

t 

V 
w 
X 

x,y,z 

(X • 

] 

Y 


Mach number 

number of half -waves in streamwise direction 
inplane loading in x-direction, positive in compression 
inplane loadii^ in y-direction, positive in compression 
number of tiles 
roots to equation (22) 
free-stream dynamic pressure 

nondimensional inplane loading parameter in x-direction, Nxb^yTr^Dj 
buckling load with no airflow 

nondimensional inplane loading parameter in y-direction, Nyb^/4^Di 

coupling-medium spring parameter defined by equations (18) 

time 

free-stream velocity 

lateral deflection 

deflection shape in x-direction 

Cartesian coordinates 

coefficients defined by equation (A5) 

compressibility factor, - 1 

plate mass per unit area 
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complex frequency coefficient, <p 

X 

dynamic -pressure parameter, 2qa^j3Di u 

Mx’^y 

Poisson^s ratio in x-direction and y-direction, respectively 

P 

free- stream air density 

<t> 

real part of complex frequency exponential coefficient 

[n] 

boundary condition matrix 

U) 

circular frequency 

c^r 

reference frequency, see equations (18) 

Subscripts: 


L 

lower plate 

le 

leading edge 

max 

maximum 

te 

trailing edge 

U 

upper plate 


Two terms, one above the other, within braces ^ ^ occur in the equations of this 
report. This notation means that two equations may be obtained by substituting first the 
upper and then the lower term. For example, equation (4b) 
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°x,U 1 

'e\ 

^ ■ '^x,U^y,u' 
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y>u ^ 2 
9y 


+ k. 
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r,te,U ax 


= 0 


at 


X = a 


ANALYSIS 

The configuration to be analyzed is shown in figure 1 and consists of an upper ortho- 
tropic plate exposed to the airflow coupled, by a linearly elastic medium, to a lower ortho- 
tropic plate. The coupling medium behaves like a Winkler foundation in that an applied 
force causes deflection only at the point of application of the force. Both plates are sim- 
ply supported along the side edges but are supported by deflectional, rotational, and tor- 
sional springs of arbitrary stiffness at the leading and trailing edges. Each plate may 
also have different material properties and inplane loadings. 

Differential Equations and Boundary Conditions 

The appropriate small-deflection equilibrium equations for the two orthotropic 
plates in the presence of inplane tensile or compressive loads may be written in the fol- 
lowing form (ref. 1); 

Upper plate 
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Lower plate 
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( 2 ) 
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The subscripts U and L refer to the upper and lower plates, respectively; for 
each plate the bending stiffnesses in the x- and y-directions are denoted by and Dy, 

respectively; Dxy denotes the plate twisting stiffness, y is the plate mass per unit 
area, and k denotes the spring stiffness of the elastic coupling medium in compression. 
The lateral loading induced by supersonic flow over the upper plate is assumed to be given 
by two-dimensional quasi -steady aerodynamic theory, so that 


2q 


dvf 




U 


8t 


( 3 ) 


where q = is the free-stream dynamic pressure and /3 = - 1 is the compres- 

sibility factor. The first term on the right side of equation (3) corresponds to a static 
loading and the second term is an aerodynamic damping term. References 5 and 6 have 
shown that for M > 1.6 and a/b from 0 to 10, use of two-dimensional aerodynamic 
theory yields flutter results in good agreement with those predicted by more exact aero- 
dynamic theory. 

According to reference 7 the boundary conditions may be written as follows: 

Upper plate 

a2w„(x,y,t) 

Wu(x,y,t) = = 0 at y = 0 and b (4a) 
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(4c) 
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Lower plate 

Replace subscript U by 


L in equations (4a) and (4b); equation (4c) becomes 
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(4d) 


The frequency-independent, linear, hysteretic formulation of structural damping 
used in reference 6 is employed in the present analysis by modifying all the terms in 
equations (1), (2), and (4) containing plate bending stiffnesses or spring stiffnesses by an 
appropriate (1 + ig) factor for simple harmonic motion. This procedure is equivalent to 
assuming orthotropic structural damping properties as well as orthotropic stiffness prop- 
erties for both plates and permits arbitrary levels of structural damping to be considered 
for each of the leading- and trailing- ec^e support springs as well as for the elastic cou- 
pling medium. 


Solution 

With introduction of the generalized aerodynamic loading terms from equation (3) 
and of structural damping and with division by the appropriate ~ 

the upper and lower plates, equations (1) and (2) become 

Upper plate 


^x,U 4 \ xy,U D * ^x,U^y,u| o 2 d 4 D 2 

8x^ \ ^1,U ’ Qx^ ^1,U 8y^ ^1,U 8x^ 


8^Wt 


8^Wt 


8^w^ 


^y,U K f \ pc 2q 

+ — LiiL }± + — X ii + Wtt - W, + + 

^1,U 8y^ ^1,U 8t^ ^1,U' ®1,U ^^1,U 


= 0 
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Lower plate 
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( 6 ) 


where the plate stiffness coefficients are defined as 
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Mx/^y 


1 - 
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(7) 


and the damping factors are defined as 

Gx = 1 + igx Gxy = 1 + igxy Gy = 1 + igy K = k(l + igi,.) 

Similarly the boundary conditions for the upper plate become 


( 8 ) 
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The boundary conditions for the lower plate become 


0^WL(x,y,t) 

WL(x,y,t) = — = 0 at y = 0 and b 
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where the complex spring coefficients are defined as 



%,le = 1 + ^^d,le 

^d,te = 1 + igd,te 



\le = 1 ^h,le 

Gf,te = 1 + iSt,te 

> 

(11) 

Gr,ie = l + igr,ie 

^r,te " 1 + igr,te 




For exponentially varying motion of the two plates the assumption of simply sup- 
ported side edges admits the followir^ lateral deflection w for the two plates as a solu- 
tion to equations (5) and (6) : 




Wy, = X,, sin — e^^ 

u u b 


V . Try i“t 
Wj = Xj sin — e^ 

L L b 


(12) 
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where X is the deflection shape in the x- direction and ^ is of the form 


? = 0 + ict) (13) 

Reference 8 has shown that for simple supports along the lateral edges there is no 
stiffness coupling between cross-stream modes. Use of two-dimensional aerodynamic 
theory precludes any aerodynamic coupling of orthogonal cross-stream modes; thus, 
modes with a single half-wave in the cross-stream direction have the lowest natural fre- 
quencies and are used to determine the dynamic stability criterion (lowest value of \ 
for which 0 > 0) for the coupled-plate system under consideration. 

Substitution of the expressions in equation (12) for w^ and w^^ into equations (5) 
and (6) results in the following two coupled fourth-order ordinary differential equations 
in and Xj^: 
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(15) 


where the primes denote derivatives with respect to x/a and where 
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2 

Similar expressions for ^L’ ®L’ "r.L’ Rjj^L, and Ry^L obtained by replacing 
subscript U with L in equations (16) and (18). 

Equations (14) and (15) may be solved simultaneously for Xy and or com- 
bined to yield a single eighth-order ordinary differential equation in Xl as follows. 
From equation (15), 


Xu - 






(19) 


Substitution of the expression in equation (19) for X^ into equation (14) yields the fol- 
lowing eighth-order linear homogenous ordinary differential equation in X^: 


Gx,U^x,l4®^ ^ ^ ^ ^ (Gx,U®L ^ ^ 






+ XAt X"' + f A„Bt + At BttIX;' + XBt X; + (BttBt - TT^SySjJX^ = 0 




( 20 ) 


A closed-form solution to equation (20) with constant coefficients may be written 


P,x/a P«x/a 


X- jA/ a. 

Xj__ = Aje + A26 


+ . . . + A0e 


Pox/a 


( 21 ) 


where the P^ satisfy the characteristic equation 


‘^x,U^x,L^i + (^x,U-^L ^x,L-^u)^f ^^x,L^i + (‘^XjU^L ^x,L®uj^f 

A T|3 .( a Ti . rt A \t»2 \ t% 


H- XA^pf + A^Bl + BuAl p2 + XB^Pi + B^B^ - J = 0 


( 22 ) 


The boundary conditions for both plates can be written in terms of the P^ and A^ 
by using equations (12) and (21) along with equations (9) and (10). When written in matrix 
form the boundary conditions become 


[n]{Ai} = 0 (23) 

Detailed expressions for the elements of the boundary condition matrix [J2 ] and their 
derivation are given in the appendix. 
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The flutter behavior of the coupled two-plate system with either aerodynamic 
or structural damping present is governed by the variation of the complex frequency 
? = ^ + io) as a function of the aerodynamic load X. Flutter of the system occurs when 
(f> becomes positive and, hence, the lowest X condition when (f> = 0 is taken to be the 
critical condition. This condition corresponds to sustained simple harmonic motion and 
represents incipient flutter for the system. 

Flutter solutions for the coupled two -plate system are obtained as follows. The 
physical characteristics of the system are specified except for the aerodynamic load X 
and the complex frequency ?. An initial aerodynamic load X is assumed and then 
incremented, and values of ? are determined for each value of X by an iteration 
scheme. The iteration on ? is continued imtil a set of roots to equation (22) is obtained 
which satisfies the boundary conditions. Since, in general, the Aj in equations (21) and 
(23) are not zero, the boundary conditions are satisfied when the determinant of the bound- 
ary condition matrix is zero. This process is continued until a value of X is reached 
which causes the real part of ? to be positive by a small value and equation (12) implies 
that the amplitude of vibration for the system increases without limit. When both aerody- 
namic and structural damping in the system are zero the real part of ? {<p) is zero for 
all values of X below some critical value. When this critical value is reached, two 
natural frequencies of the system coalesce into complex conjugates, so that positive val- 
ues of the real part of ? (<(>) occur and equation (12) again implies motion with growing 
amplitude, that is, flutter. 

To obtain flutter boundaries for the coupled two-plate system a computer algorithm 
was written to obtain values of ? as a function of X. The algorithm is written in 
FORTRAN IV for the CDC 6000 series computers at the Langley Research Center and 
makes use of library subroutines to numerically extract the roots of the eighth-order 
complex polynomial (eq. (22)) and to evaluate the determinant of the complex boundary 
condition matrix. In the algorithm, values of ^ which satisfy equation (23) are deter- 
mined for given values of X by using a Newton-Raphson iteration scheme wherein a new 
value of ? is determined from the current value by the following expression: 

|n(?i)|AC 

iK'i)] - * *?) 

where is a small increment in ? used to numerically calculate the first derivative 
of the determinant of the boundary condition matrix with respect to 
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RESULTS AND DISCUSSION 


Because of the large number of variables involved, no attempt was made to conduct 
extensive parametric studies of the flutter behavior of coupled two-plate systems. How- 
ever, results are presented to show that when aerodynamic and structural damping effects 
are ignored and the two plates are identical with respect to geometry, inplane loadii^, 
ec^e conditions, and material properties, flutter boundaries can be obtained in terms of 
the general parameters A and B as was done in references 8 and 9 for single plates. 
Some effects of structural damping on the flutter of coupled identical plates for a/b = 4 
are also presented. In addition, results from the current analysis are compared with 
results from the modal solutions of references 1 and 2 to assess the validity of these 
solutions. Two thermal protection systems proposed for space shuttle which can be 
modeled as coupled two-plate systems are analyzed to determine their flutter character- 
istics. For each system the stiffnesses of the upper and lower plates differ by 2 orders 
of magnitude. The plate with low stiffness is exposed to the flow in one system and in the 
other the plate with high stiffness is exposed to the flow. A discussion of the relative 
flutter behavior of the two systems is also presented. Finally, the influence on flutter of 
single plates of unequal support springs at the leadii^ and trailing plate edges is discussed 
briefly. 


Flutter Boundary for Identical Coupled Plates 

When aerodynamic and structural damping are neglected and when the two plates 
are identical with respect to geometry, inplane loading, and material properties, equa- 
tions (16) and (17) indicate that the respective A and B parameters for the two plates 
are equal. For these conditions it is possible to obtain general flutter boundaries for the 
coupled system for specific stiffness values of the elastic coupling medium and specified 
boundary conditions. Such flutter boundaries are then applicable to a wide variety of 
panel configurations, material properties, and inplane loadings. 

Although it is possible to generate solutions for various values of the spring stiff- 
ness of the coupling medium and for various boundary conditions, only one such solution 
with both plates simply supported and a coupling- medium spring parameter of 
Su = Sl = 10 (representative of a lightly coupled system) was obtained to illustrate the 
nature of such solutions. Figure 2 shows a comparison of the resulting flutter boundary 
for a coupled two-plate system with a similar boundary for a single plate. The boundaries 
are shown in terms of as a function of A. The anomalous X = 0 flutter points 

occur whenever A has a value that causes the two vibration modes which coalesce flut- 
ter to have equal natural frequencies for no airflow. For the coupled-plate system the 
larger number of X = 0 flutter points results from the fact that the coupled system has 
twice the number of natural vibration frequencies associated with a single plate. 
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Eight natural frequencies of the coupled system are shown in figure 3 in which the 
frequency parameter B is shown as a function of A. The solid curves correspond to 
vibration modes of the two plates which are in phase, the dashed curves correspond to 
vibration modes of the two plates which are 180® out of phase, and m denotes the num- 
ber of half-waves in the stream direction. Shaded regions between any two curves indi- 
cate that those modes contiguous to the shaded regions on the plot coalesce to give the 
lowest value of X for flutter. As usual in panel flutter, adjacent odd and even modes 
coalesce for flutter. However, since the coupled-plate system has both in-phase and out- 
of-phase vibration modes, many possibilities exist for coalescence of multiple pairs of 
modes which have nearly the same degree of frequency separation, and it is not obvious 
which pair will give the lowest flutter value of X. Comparison of figures 2 and 3 indicates 
that the X = 0 flutter points coincide with crossings of frequencies which coalesce for 
flutter and that the multiplicity of these points results from the double set of frequencies. 

Although variation of the coupling-medium spring parameter S is not expected to 
change the basic character of the flutter boundary for the two-plate system, such changes 
will influence the separation of the in-phase and out-of-phase natural frequencies for the 
system and hence the A values at which the X = 0 flutter points occur. Reference 6 
shows that inclusion of aerodynamic and structural damping in the theory can remove 
the X = 0 flutter points for single plates. Calculations were made to determine whether 
the same effects hold for a coupled two-plate system. However, since inclusion of damp- 
ing precludes a flutter solution in terms of the general parameters A and B, it is nec- 
essary to pick a specific plate configuration. Thus, an a/b = 4 coupled-plate system is 
used to illustrate the effects of aerodynamic and structural damping on the flutter of cou- 
pled plates. Flutter boundaries for such a system with = Sl ” shown in fig- 

ure 4 in which X^/^ is shown as a function of the inplane loading parameter R^. The 
dashed boundary is for no damping and the solid boundary is for structural damping coef- 
ficients gx, gxyj Sy 0-01 in the plates and aerodynamic damping corresponding 
to sea level conditions. The aerodynamic damping is based on a 0.051-cm (0.02-in.) 
thick aluminum plate and sea level air properties and was calculated from the following 
equation from reference 6: 




The damping has a strong smoothing effect on the flutter boundary and completely removes 
the X = 0 flutter points. Although on the basis of these a/b = 4 results, it appears that 
damping has a large stabilizing effect and removes the X = 0 flutter points for both one- 
and two-plate systems, such is not always the case. For example, reference 10 shows 
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that structural damping is destabilizii^ for built-up plates with low a/b ratios (approx- 
imately 1). Limited results concerning such behavior are discussed in a subsequent 
section. 


Comparison of Modal and Closed- Form Results 

References 1 and 2 present results from two-mode GalerMn solutions for a coupled 
two -plate system with a/b = 1, simple support boundary conditions, and a value of the 
coupling -medium spring parameter Sy = Sl = S = 20. Reference 2 also presents a flut- 
ter boundary for S = 10. In figure 5(a), results from both analyses are compared with 
closed-form results from the present analysis. Flutter boundaries for identical plates 
with identical inplane loading are shown in terms of A as a function of the inplane load- 
ing parameter Rx. These results are for S = 20. Closed-form results are shown by 
the solid curve and modal results are shown by the dashed curve (ref. 1) and by the cir- 
cular symbols (ref. 2). A closed-form boundary for a single plate is shown for reference. 
As mentioned in reference 2, for some as yet unexplained reason the two modal solutions 
from references 1 and 2 do not agree, although both employ the same representation of 
the mode shapes. The two -mode GalerMn solutions have the correct trends but tend to 
be unconservative for some values of Rx- Results for S = 10 are shown in figure 5(b) 
and again the two-mode solution has the correct trend, but it shows a tendency toward 
conservatism in some regions particularly as Rx increases negatively. The results of 
figure 5 indicate that as S increases from 10 to 20, the first A = 0 flutter point shifts 
to a negative (tensile) value of Rx. As mentioned in reference 2, inclusion of aerody- 
namic damping will remove the anomalous A = 0 flutter point; but as indicated on fig- 
ure 5(a), the boundary continues to exhibit a large dip in the vicinity of Rx = -7. Calcu- 
lations from the closed-form solution showed the same trends, and when structural 
damping was included in the calculations, no appreciable effect was found at the A = 0 
flutter points; however, the structural damping was found to be destabilizing for other 
values of Rx- As mentioned previously, this destabilizii^ effect has been noted for sin- 
gle plates in reference 10 and was not pursued further in the current investigation. 

Reference 1 also presents results for two geometrically identical plates with only 
one plate subjected to inplane loading. These results are compared with closed-form 
results in figure 6. Again a/b = 1, S = 20, and both plates are simply supported. Fig- 
ure 6(a) shows flutter boundaries for the upper plate only subjected to inplane loading and 
figure 6(b) shows similar boundaries for the lower plate subjected to inplane loading. 

Solid curves denote closed-form results and dashed curves denote two-mode results. In 
this instance the two- mode solutions are for the most part unconservative and in the 
region of negative Rx do not exhibit correct trends of the closed-form solution. Exam- 
ination of figure 7, which shows the natural frequencies of the coupled two-plate system 
with only one plate subjected to inplane loadii^, reveals the reason for this behavior. 
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Shaded regions denote that those modes contiguous to the shaded regions coalesce for 
flutter. Since the two-mode solution is restricted to the first two modes, it can give 
reliable results only when these modes coalesce for flutter. Thus, the disparity between 
the exact and two-mode solutions results from the fact that the boundary for the upper 
plate subjected to inplane loading (fig. 6(a)), except for < -19 and Rx > 6.5 (fig. 7), 
is governed by a coalescence involving the third natural frequency. When the lower plate 
is subjected to inplane loading (fig. 6(b)), the flutter boundary from Rx = -13 to 0 is 
governed by a coalescence involving modes 1 and 2, and although the two -mode solution 
is unconservative, the trend is correct. Outside the region Rx = -13 to 0, mode 3 is 
again involved and the disparity between the two-mode and closed-form solutions worsens. 

Application to Proposed Thermal Protection Systems 

Flexible mat system .- In reference 11 a thermal protection system (TPS) is pro- 
posed which consists of a flexible mat fabricated from fibrous ceramic felt layers covered 
by a dense ceramic coating bonded to a load carrying primary structure. Since no flutter 
results for this TPS were presented in reference 11, the present study was undertaken. 

The TPS was modeled as a system of springs connecting the surface covering to the sub- 
structure, and calculations were made to determine the flutter characteristics for such 
TPS. The system, shown in figure 8, is similar to that considered in figure 6(b) with the 
exception that the plate exposed to the flow has a much lower flexural stiffness than the 
unexposed plate and the ends of the exposed plate are not necessarily simply supported. 

Flutter boundaries which were calculated for this system on the assumption that it 
behaves as two spring-connected beams are also shown in figure 8. The solid curves 
represent flutter boundaries for beams of finite stiffness. The upper boundary results 
from the assumption that the leading and trailing edges of both beams are simply sup- 
ported. The lower boundary corresponds to the assumption that the leading and trailing 
edges of the upper beam are supported by deflectional springs with a spring stiffness 
equal to that of the coupling medium. The two horizontal dashed curves correspond to 
similar boundary conditions but with the lower beam of infinite bending stiffness. 

Figure 9 shows the natural frequencies of the two-beam system as a function of 
the inplane loading in the lower beam. The solid curves are for simple supports and in- 
phase motion, the dashed curve is the first out-of-phase mode for simple supports, and 
the dot- dashed curves are for spring supported edges and in-phase motion. Shaded 
regions between any two curves indicate that those modes contiguous to the shaded 
regions coalesce for flutter. Again the difference in the two flutter boundaries is related 
to the modes which coalesce for flutter. Along the boundary for simple supports, in-phase 
modes 1 and 2 coalesce for flutter up to Rx/^x,cr = 6.5; between Rx/Rx,cr = 6.5 and 
approximately 0.67, in-phase modes 2 and 3 coalesce; and beyond Rx/Rx,cr “ 6.67, in- 
phase mode 2 and out-of -phase mode 1 coalesce for flutter. For the boundary cor re- 
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sponding to spring supports, in-phase modes 2 and 3 coalesce along the entire boundary 
to give the flutter condition. For this structure the upper beam has low flexural stiffness 
and the lower beam has high flexural stiffness (2 orders of magnitude greater than the 
upper). For both sets of boundary conditions it is interesting to note that if the stiffness 
of the lower beam is considered to be infinite, conservative predictions of X result 
which are in good agreement with the corresponding minimum values of X obtained for 
a lower beam of finite stiffness. This suggests that for configurations in which the lower 
beam carries the inplane loading and is very stiff in relation to the upper beam, use of the 
simpler solution may be acceptable. 

Reusable su rface iji^lation.- Reference 3 presents a modal solution for the flutter 
behavior of a proposed space shuttle TPS consisting of relatively thick ceramic tiles 
mounted on a soft viscoelastic foundation which is bonded to a primary load carrying 
metallic structure. The concept is usually referred to as a reusable surface insulation 
(RSI) TPS. A Rayleigh-Ritz flutter solution employing 4 free-free beam streamwise 
modes and 1 free-free beam cross-stream mode per tile, and 12 simply supported plate 
streamwise modes and up to 3 simply supported plate cross-stream modes for the metal- 
lic primary structure was obtained. 

The structural model in reference 3 for the viscoelastic foundation neglects viscous 
effects of the foundation material but is classified as a shear model (in contrast to the 
Winkler model) in that a point load gives rise not only to a deflection at that point but to 
surrounding points as well with an exponential decay of deflection away from the point of 
application of the force. In reference 3, the Winkler model and shear model were found 
to give essentially the same results for a single tile and for streamwise multiple tile 
arrays when the viscoelastic foundation was assumed to be cut between tiles. The cut in 
the foundation is a part of the space shuttle design; hence, the difference in foundation 
models does not preclude valid comparisons of results from the present analysis with 
results from reference 3. Other differences also exist between the two analyses. For 
instance, the analysis of reference 3 is applicable to multitile arrays and assumes the 
tiles to be free on all edges. The current analysis, however, is restricted to a single 
two-plate combination which is simply supported along the side edges. Thus, the two 
analyses can be compared only on a limited basis. One valid comparison is for a/b = 0. 
Reference 3 gives a flutter value for X of 65 for a two-plate system where the RSI tile 
has a flexural stiffness of 52 N-m (458 lb-in.), the foundation has coupling -medium spring 
parameters of Sy = 32 and Sl = 1971, the lower plate has a flexural stiffness of 
0.75 N-m (6.67 lb-in.), the plates are 51 cm (20 in.) long, and sea level aerodynamic 
damping is included. The present analysis gives a value of X of 56 for similar 
conditions. 

An additional comparison with results from reference 3 is shown in figure 10 in 
which flutter boundaries are given for an a/b = 4 coupled-plate system. Since refer- 
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ence 3 assumes the tile to be free on all edges, the upper plate was modeled for the cur- 
rent analysis as having stiffness in the streamwise direction only and hence is a beam- 
type representation of the actual tile. The flutter boundaries shown in figure 10 are pre- 
sented in terms of X as a function of Rx/l^x,cr ^x,cr value of the 

inplane loading parameter required to buckle the lower plate alone. Structural details of 
the coupled system are shown in the inset in figure 10. Two boundaries from reference 3 
are shown; the boundary labeled n = 1 is for a single tile, and the curve labeled n = 4 
is for a streamwise array of four tiles. The boundary labeled closed-form solution is 
from the present analysis and agrees best with the n = 4 boundary from reference 3. 

The apparent lack of agreement between the closed-form and modal results for n = 1 
may be due to the relatively small number of modes used to describe the motion of the 
tile. Two rigid-body and two bending modes may not be sufficient to adequately describe 
the motion of the tiles for n = 1. 

Results from the closed-form solution indicate that, although the mass of the upper 
plate is approximately twice that of the lower plate and the stiffness of the upper plate is 
two orders of magnitude greater than that of the lower plate, the coupling between the two 
plates is sufficient to cause the upper plate to respond predominately in bending rather 
than as a rigid body. The better agreement between the closed-form n = 1 and modal 
n = 4 flutter boundaries can be attributed to the closeness of the mode shapes for the 
respective solutions. Figure 11 shows mode shapes for the first two modes of the cou- 
pled system. Since those are the modes which coalesce for flutter, it can be expected 
that the two solutions might give similar flutter results. Reference 3 points out that the 
salient feature of this coupled system is that although the upper plate has a greater mass 
and a much higher stiffness than the lower plate, it is the flexibility of the lower plate 
which governs the response of the system. In fact, the addition of the upper plate and 
elastic coupling medium actually results in an increase in the dynamic pressure required 
to flutter the coupled system over that required to flutter the lower plate alone 

^q//3coupled ^/^single ^ psi)^. The opposite behavior 

occurs for the situation shown in figure 8 in which the upper plate has a much lower flex- 
ural stiffness (two orders of magnitude) than the lower plate ^q/^coupled " 

(121 psi) and q//3singie " MPa (15 290 psi)j. In addition, the presence of the RSI 

tiles and viscoelastic foundation increases the buckling load for the coupled system over 
that for the lower plate alone. 

Single plates with unequal supports .- As mentioned in reference 4, metallic thermal 
protection systems have been proposed which can have rows of continuous supports with 
widely varying flexibilities normal to the plane of the panel. To illustrate the effects on 
flutter of such unequal deflectional supports, calculations were made using the present 
analysis for single plates with large streamwise bending stiffness and essentially zero 
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cross-stream bending and twisting stiffnesses. Figure 12 shows stability boundaries in 
terms of X as a function of spring stiffness ratios k^j ^g/k^j from 0 to 1.0 and 
kd te/^d le stability boundaries are for constant values of the spring 

stiffness parameter dashed curve, shown for comparison, 

is the value of X for infinite spring stiffness or simple supports. 

The curves show large reductions in X from the value associated with simple sup- 
ports as the edge spring stiffness is decreased. The effect of unequal supports appears 
to be a function of ^g. For example, when ^g is large (e.g., 100), the effect of 
reducing the trailing-edge spring stiffness is small as shown for k^j j-g/k^j^^g < 1. How- 
ever, when small (e.g., 1), a large reduction in the flutter X occurs as the 

trailing-edge spring stiffness decreases. Additionally, when ^g is small, large val- 
ues of the trailing-edge spring stiffness ^kjj^^g/lqj^te l) introduce divergence (loss 
of static stability) which also results in large reductions in X. 

On the basis of the results presented in figure 12, it appears that in the design of a 
TPS which employs flexible edge supports to alleviate thermal stresses, adequate atten- 
tion must be given to the design of the supports to preclude the drastic reductions of aero- 
elastic stability margins associated with unequal supports. 

CONCLUDING REMARKS 

Flutter of two parallel elastically coupled plates was investigated analytically. A 
closed-form solution including both aerodynamic and hysteretic structural damping was 
presented for supersonic flutter of flat orthotropic plates coupled by an elastic medium. 
Both plates were simply supported along the side edges but were supported by deflectional, 
rotational, and torsional springs of arbitrary stiffness at the leading and trailing edges. 
Two-dimensional quasi-steady aerodynamics was used in the analysis. 

Because of the large number of variables involved, no attempt was made to conduct 
extensive parametric studies of the flutter behavior for coupled two-plate systems. How- 
ever, selected results were presented to illustrate the general flutter nature of such sys- 
tems. Similar to results for sir^le plates, anomalous zero-dynamic-pressure flutter 
points were predicted for the two -plate system whenever the frequencies which coalesce 
for flutter were equal for no airflow. Inclusion of aerodynamic and structural damping in 
the calculation removed the computed anomalous zero-dynamic-pressure flutter points for 
the two-plate system for high length-width ratios (approximately 4) but had little effect for 
low lei^th-width ratios (approximately 1). 

A comparison of flutter results was made for two thermal protection systems pro- 
posed for space shuttle. The two systems were similar in that for each the coupled plates 
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differed in flexural stiffness by two orders of magnitude. The two systems differed in 
that for one the plate of higher flexural stiffness was exposed to the airflow and in the 
other the plate of lower flexural stiffness was exposed to the airflow. The comparison of 
flutter results revealed that for coupled two-plate systems with different flexural stiff- 
nesses in the two plates, maximum flutter margins occur when the plate with the higher 
flexural stiffness is exposed to the flow. 

Finally, on the basis of flutter results presented for a sii^le plate with unequal 
deflectional supports, it appears that for thermal protection systems which employ flex- 
ible edge supports to alleviate thermal stresses, adequate attention must be given to the 
design of the supports to preclude the drastic reductions in aeroelastic stability margins 
associated with unequal supports. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
November 4, 1975 
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APPENDIX 


BOUNDAHY CONDITION MATRIX 

The boundary conditions are given by equations (9) and (10) in the main text. Equa- 
tion (12), which is assumed to satisJ^ equations (5) and (6), satisfies exactly boundary con- 
ditions (9a) and (10a) and is repeated for convenience; that is, 


w-Q = sin ^ e^^ 


Wt - Xt sin — 

L L b 

These expressions modi:fy the remaining boundary conditions (eqs. (9b), (9c), (10b), and 
(10c)) as follows; 

Upper plate 


,r2a2 





at X = 


-G. 


d,Ze,U^d,Ze,UU^ \ 


Xt 


'+‘^d,te,U^d,te,Uj 


"U ""LA K U 

■t,te,U^t,te,Uj 



^x,U^U > 


Trar/or- 


^-sr TtAtt rr “ TT ^T1 


= 0 at X = < 
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where 


^d,te,U 

^r,Ze,U 

^t,Ze,U 


^d,Ze,U^^ 

^1,U 

^r,^e,U^ 

^1,U 

^,^e,U^ 

°1,U 


^d,te,U 

^r,te,U 

*^,te,U 


^d,te,U^^ 

°1,U 

^r,te,U^ 

°1,U 

^,te,U^ 

^1,U 


These expressions for the lower plate may be obtained by replacing subscript U with L. 
Finally , 


p _ ^d,te,U^^ 

^d,ie,U - Dj L ' 


:p _ l^d,te,U^^ 

^d,te,U - Di L ' 


Equations (19) and (21) from the main text (repeated here for convenience) are 


Xtt = (G„ T + At xV + Bt X 


U “ “1“ rx,L^L ^ 


Pi x/a Po x/a Pq x/a 

Xj^ = Aje ^ + A2e ^ + . . . + Age ° 

With these equations the boundary conditions given by equations (Al) and (A2) may be 
written as follows: 

Upper plate 



p2 _ n (]I^ 
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where 




For j = 1, 2, 3, . . .,8. The eight equations given by equations (A3) and (A4) can be 
combined and written in matrix form as follows: 


[ S2]{A} = 0 


(A3b) 


(A4) 


(A5) 


(A6) 
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where [S2] isa8x8 square matrix and is a 8 x 1 column matrix. The elements 

of [^2] are given by the following expressions: 
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for j = 1, 2, 3, . . ., 8. 
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Figure 1.- Coupled plates, coordinate system, and spring support system. 



Figure 2.- Flutter boundaries for parallel elastically coupled identical plates and for 
a single plate; S = 10; no structural and aerodynamic damping. 
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Figure 4.- Flutter boundaries for simply supported identical plates, a/b = 4; S = 10 
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(a) S = 20. 

Figure 5.- Closed- form and two-mode flutter boundaries for simply supported 

identical plates, a/b = 1. 



Flow 



(b) S = 10. 

Figure 5.- Concluded. 





(a) Upper plate subjected to inplane loading. 

Figure 6.- Closed- form and two-mode flutter boundaries with one plate subjected to 

inplane loading. S = 20. 
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Figure 7.- Natural frequencies for elastically coupled, simply supported plates 
with one plate subjected to inplane loading, a/b =1; S = 20. 
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Figure 8.- Flutter boundaries for proposed TPS. 





Figure 9.- Natural frequencies for proposed TPS. 
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